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Levy-type diffusion on one-dimensional directed Cantor Graphs 
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Levy-type walks with correlated jumps, induced by the topology of the medium, are studied on 
a class of one-dimensional deterministic graphs built from generalized Cantor and Smith- Volterra- 
Cantor sets. The particle performs a standard random walk on the sets but is also allowed to 
move ballistically throughout the empty regions. Using scaling relations and the mapping onto the 
electric network problem, we obtain the exact values of the scaling exponents for the asymptotic 
return probability, the resistivity and the mean square displacement as a function of the topological 
parameters of the sets. Interestingly, the systems undergoes a transition from superdiffusive to 
diffusive behavior as a function of the filling of the fractal. The deterministic topology also allows 
us to discuss the importance of the choice of the initial condition. In particular, we demonstrate 
that local and average measurements can display different asymptotic behavior. The analytic results 
are compared with the numerical solution of the master equation of the process. 
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I. INTRODUCTION 



Stochastic processes characterized by Levy walks are 
relevant in many physical phenomena, ranging from con- 
densed matter, biological systems and ecology to trans- 
port in turbulent fluids and in porous media. In recent 
years there has been a growing interest in the field (see 
Ref. [l| for a recent overview). However, direct compar- 
ison of theoretical predictions with experimental data in 
fully controlled experimental conditions is often unfeasi- 
ble. As a matter of fact, evidence of anomalous transport 
and diffusion is often indirect and experimental setups 
allowing a direct and tunable study of such phenomena 
are hardly available. A breakthrough in the field was 
achieved in a recent experiment, by means of an engi- 
neered material where light rays were demonstrated to 
perform a Levy walk-type of diffusion Q . This new class 
of materials is built by packing glass microspheres with 
diameters following a Levy distribution, and then filling 
the space between the spheres with strongly scattering 
nanoparticles. As light is not scatterered within the mi- 
crospheres, anomalous diffusion arises as a consequence 
of the power-law distribution of their diameters. The re- 
sults of the experiments were indeed interpreted using a 
model of Levy walks [3, H| for photons, and the new ma- 
terial was named Levy glass. An important feature of the 
experimental samples is that the walk is correlated, and 
the correlation is induced by the topology of the quenched 
medium. Light that has just crossed a large glass micro- 
sphere without being scattered has a high probability of 
being backscattered at the subsequent step undergoing a 
jump of similar size. Studies on models of Levy flights 
have evidenced that spatial correlations in jump prob- 
abilities have a deep influence on the diffusion proper- 
ties 043- Here, however, we consider Levy walk models, 



where particles either move ballistically or are scattered 
by impurities, and do not perform long range instanta- 
neous flights. Recently, one-dimensional models for Levy 
walks, where the correlation is induced by an underly- 
ing quenched and correlated random environment, have 
been studied. One-dimensional models represent simpli- 
fied systems and may not compare quantitatively with 
real experiments. Nevertheless they allow for an exact 
solution of the dynamics. The recent studies focused, 
respectively, on the mean square displacement in a Levy- 
Lorentz gas @ , and on the conductivity and transmission 
through a chain of barriers with Levy-distributed spac- 
ings [9j. Interestingly, in simple random walks these dif- 
ferent aspects are not independent and they can be con- 
nected by assuming standard scaling relations for the ran- 
dom walk probabilities [10(, and exploiting the analogy 
between the associated master equation and the Kirch- 
hoff equations 11| . When an analogous scaling holds for 



Levy walks, then it can be applied to relate the scal- 
ing exponents of the return probability, the mean square 
displacement and the effective resistance of the samples. 
Another subtle point in random models is represented 
by the choice of the starting site, as discussed in Ref. [|| , 
and it would be interesting to have a direct control on 
this problem. 

In this paper, we will present a new class of one- 
dimensional models for correlated Levy walks, that will 
allow us to investigate in detail the scaling and the re- 
lation between transmission and mean square displace- 
ment, as well as the dependence on the starting point 
through average and local diffusion properties. More 
specifically, we will consider diffusion across deterministic 
fractals, namely generalized Cantor and Cantor-Smith- 
Volterra sets (the latter being the simplest examples of 
fat fractals). The walker dynamics is such that it per- 
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forms a standard random walk step on sites belonging to 
the set, while moving ballistically (without changing di- 
rection) otherwise. In the following we will refer to these 
two types of sites as "bidirectional" or "unidirectional" , 
respectively. Levy-type diffusion may thus arise as the 
longer steps will be distributed according to the voids of 
the Cantor set, much in the same spirit of the Levy glass 
experiment [2|. 

The model can be mapped onto a random walk on a 
suitable directed graph and it allows to obtain analyti- 
cal results by extending the usual scaling arguments. In 
particular, the scaling of the resistivity with the length 
of the graph yields the exponents for the asymptotics 
of the return probability and the mean square displace- 
ment. Interestingly, the system undergoes a transition 
from superdiffusive to diffusive behavior as a function of 
the density of the fractal. A similar transition has already 
been evidenced for Levy flights in presence of quenched 
noise @ and also in the model of Ref. @ . 

A subtle effect of quenched disorder is the dependence 
of the observables on the choice of the walker starting site 
Q. An insight on this issue is clearly relevant for guid- 
ing experimental investigations. Due to the simple deter- 
ministic topology, our model allows for a clear analysis of 
such issue. In particular, the choice of a random starting 
point corresponds to a measure of quantities averaged 
over all possible starting sites. On inhomogeneous struc- 
tures, such averages can display a different asymptotic 
behavior [l2l ] with respect the corresponding local quan- 
tities, obtained by taking into account a specific starting 
point. This is indeed the case for generalized Cantor 
sets: scaling relations hold for local quantities but are 
violated for average ones. In particular, in the local case, 
the random walker shifts from a superdiffusive to a dif- 
fusive motion, as a function of topological parameters 
tuning the void density. In the average case the behav- 
ior is even richer: a change of the void density and size 
induces a transition among three different regimes char- 
acterized by ballistic, superdiffusive and standard diffu- 
sive motion, respectively. In this case, which appears to 
be the closest to the experimental settings, adopting a 
method similar to Ref. [8(, we give a rigorous bound on 
the mean square displacement and discuss the violation 
of scaling relations. 

The paper is organized as follows: In Sec. [n] we intro- 
duce the topology and the dynamics of the model. In 
Sec. IIIII we first review the scaling relations connecting 
random walk properties and the resistivity and derive the 
exact results for the resistivity and for the local asymp- 
totic behavior of the walker itself. We also check, via the 
numerical solution of the master equation, the scaling hy- 
pothesis for local quantities. In Sec IIVI we turn to the 
average case and provide a rigorous lower bound for the 
exponent of the mean square displacement, as a function 
of the distribution of the voids. We show under which 
condition the bound is fulfilled exactly. Also in this case 
we compare our analytical results with the numerical so- 
lution of the master equation. Some concluding remarks 
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FIG. 1: (Color online) Upper figure the discrete Cantor set 
in grey (fuchsia) the sites which have been removed (voids). 
Lower figure the Cantor graph associated to the Cantor set 
at generation G = 3 where undirectional bubbles have been 
introduced to describe ballistic motion in the voids. The site 
i = 1 is the origin of the graph. 

are given in Sec. [V] 



II. A RANDOM WALKER ON GENERALIZED 
CANTOR GRAPHS 

A. Random walk across Cantor Sets 

The standard Cantor set is built by removing the mid- 
dle thirds of a linear chain and iterating the rule in every 
segment. The resulting set has zero measure, with fractal 
(Hausdorff) dimension = j^jf- The procedure can 

be easily generalized: let us denote by n s the number of 
parts in which we split the segment (n s odd). Remov- 
ing parts in even positions and iterating the process, the 
measure of the remaining set is still zero and its fractal 
dimension is 



If one builds the Cantor set by growing the fractal, then 
in this notation n r = " , ' 2 t " 1 is the number of replicas of the 
generation G — 1 in the generation G. The integer n r will 
be used in the following to identify a given set. One can 
introduce the discrete version of the above procedure, as 
in the upper part of Fig. [T] where the dark sites represent 
the Cantor set and the grey (fuchsia) sites have been 
removed and correspond to voids. 

Let us now introduce a random walker on such fam- 
ily of Cantor set. The walker can reside on any site of 
the lattice. If it is on the void regions, it moves bal- 
listically without changing direction: when the random 
walker enters a void from left (right), it keeps going left 
(right) without being scattered until it reaches the oppo- 
site edge of the void itself. If instead the walker is on the 
sites belonging to the Cantor set, it performs a standard 
random walk with 1/2 probability to move right or left 
at the next step. 

Denoting with P^(t) the probability that a walker 
started at site i arrives at time t at site j, respectively 
with positive or negative velocity (i = 1, . . . , N tot can be 
both a site of the Cantor set or a void site, see upper part 





FIG. 2: Three examples of Cantor graphs of generation 2 with 
n r = 2,3,4. 



FIG. 3: The generation G = 3 of the graph with n u = 3, 4, 5, 
in the n r = 2 case. 



of Fig. [T]), then the master equation of the corresponding 
process reads 

P+(t + 1) = Tj-xP+^t) + (1 - Tj_x)P^_ x (t), 
Pi](t + 1) = (1 - T i+1 )P+ +1 (t) + Tj+xP^t). (2) 

The Cantor structure is therefore defined by the trans- 
mission coefficients 



1/2, if j is on the Cantor set 



1. 



on the voids 



(3) 



The above choice for Tj yields indeed a ballistic motions 
on the voids and a random-walk behavior on the Cantor 
set. The choice Tj = T is readily recognized to corre- 
spond to the standard persistent random walk [l3|. In 
this respect, our model can also be regarded as a random 
walk with a site-dependent persistence (see e.g. [13] and 
references therein). 



B. Cantor Graphs from Cantor Sets 

Let us introduce an equivalent approach, based on 
graph theory. We associate to the Cantor set a directed 
Cantor graph, obtained by putting undirected links in 
every "solid" segment and replacing the voids with bub- 
bles of appropriate length whose sites are connected by 
directed links (see Fig. Q] for an explicit construction). 
Some examples of graphs at different n r are shown in 
Fig. [5] In this framework, the model consists of a ran- 
dom walk on such a directed graph. The adjacency ma- 
trix Aij of the graph is defined as Aij = Aji = 1 if sites 
i and j are connected by an undirected link, = 1 if 
the site i is connected to site j by a link directed from i 
to j and otherwise. The outgoing coordination number 
Zi = Aij is one in the voids and two in the solid seg- 
ments. The transition probability for the random walk 
is therefore: = Aij/zi. 

The structure of the voids in the Cantor graph strongly 
influences the motion of the random walker, so let us 
now introduce a further generalization. In the previous 
examples, the lengths Lfc of the unidirectional segments 
are Lk = n k s , k = 1, . . . , G — 1, where G is the generation 
of the graph. However, we can choose an arbitrary n u 



(n u integer > 2), so that the length of the unidirectional 
segments becomes 



Li 



k = 1, 



■,G- 1, 



(4) 



as is shown in Fig. [3] 

From Eq. (QJ, we can compute Nt,(G,n r ) and 
N u (G,n r ,n u ), the number of bidirectional and unidirec- 
tional sites respectively (those plotted in black and fuch- 
sia in Fig. [IJ in a graph of generation G: 



if n u ^ n r 



+ (n r — l)Gn^ 1 + 1 if n u = n r 



(5) 



(6) 



is 



Accordingly, the total number of sites is Not {G, n r ,n s ) 

(2n r -n„-l)ng-(n r -l)ng 1 jf ^ 

n r —n u a / i 

V/„/ = [ (7) 
+ (n r - l)G?i^" 1 + 1 iin u =n r 

From Eqs. ([5]), © and J?]), it follows that the measure 
of both sets Nf, and N u is strictly positive if n u < n r . 
This choice of the parameters corresponds to the so called 
fat fractals. An example of this kind of structure is 
the Smith- Volterra-Cantor set. Conversely, we will name 
slim fractals the graphs with n u > n r , where the measure 
of the solid segments is vanishing. 

It is important to notice that, even if the Cantor set has 
a fractal dimension di se *^ < 1, Eq. ((T|), the directed Can- 
tor graph has always a fractal dimension df = 1. Indeed, 
since the voids have been filled with one-dimensional bub- 
bles, the number of sites within a distance r from a given 
site grows as ~ r df , with df = 1. 



III. DIFFUSION ON CANTOR GRAPHS: 
LOCAL BEHAVIOR 

A. Asymptotic properties of random-walks, 
spectral dimensions and Einstein relations 

On the directed Cantor graph, we can now define the 
standard random walks quantities. Their asymptotic 
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behavior are expected to be described by power laws, 
with characteristic exponents, usually related only to the 
large-scale topology of the structure [l2[ . 

Let Pij(t) be the probability that a walker started at 
site i arrives at j in t steps, P%j{t) — P^(t) + Pjjit). 
In particular, Pu(t) is the random walk autocorrelation 
function i.e. the return probability to the startin g p oint 
after t steps. At large times t it is expected that |15| 



Pu(t)~r 



(8) 



where d s often coincides with the spectral dimension, 
which also rules the low eigenvalues region of the spec- 
trum for the discrete Laplacian on the graph [lfij]. In 
the case at hand, being the graph directed, d s simply 
describes the asymptotic behavior of the return proba- 
bility, but we will refer to it as a "generalized" spectral 
dimension as well. 

The mean square displacement from the starting site i 
after t steps is: 




(9) 



where the distance between the sites i and j. If the 
starting site corresponds to the origin of the Cantor graph 
(i = 1 in Fig. [1]) and periodic boundary condition are 
chosen, as in all our numerical solutions, then (x^ = 
for simmetry reasons. The exponent 7 defined by the 
asymptotic behavior 



(10) 



classifies the diffusive properties of the random walker: 
7 = 1 corresponds to usual diffusion, 7 > 1 to superdif- 
fusion, and 7 < 1 is typical of a subdiffusion. The value 
7 = 2 characterizes a ballistic motion where the random 
walker displacement grows linearly with time. 

It is well known that on regular structures, the dy- 
namic exponents ruling the return probability and the 
mean square displacement are not independent, i.e. for 
systems where scaling holds [Ioj |. 



7 = d s /df 



(11) 



In the Appendix we will show that such a relation still 
holds under a generalized scaling hypothesis that takes 
into account the underlying fractal structure of the Can- 
tor graph. More precisely, we assume that Pij (t) depends 



on j only through the distance Xi 
i.e. Pij(t) = Pi(r,t) and that 



between i and j , 



P;(r, t) 



(12) 



where £(t) is the correlation length of the system and g(-) 
is a periodic function of a logarithm in arbitrary base x. 
The master function fa can depend on the initial site. 



The scaling factor t~ ds / 2 ensures consistency with the 
definition of the return probability, Eq. ©, as seen by 
setting j = i (r = 0) in Eq. (Tl2|) . Moreover, log-periodic 
oscillations are superimposed to the leading exponential 
behavior as it is expected on fractals [2(J. In the Ap- 
pendix we also show that, up to log-periodic corrections, 
(x?) ~ £ 2 (t) and that 



£(t) = t^s 



(13) 



In general Pi (r, t) presents the above scaling form for 
large enough times. In our calculations we verified that if 
the walker starts from the origin of the structure, scaling 
is realized in a few steps. A different choice can give rise 
to very long transients where Pijit) may also depend on 
the direction and not only on the distance r alone. The 
net effect is a nonvanishing drift (xi(t)) ^ for short 
times. 

A simple analogy between the master e quat ion of the 
random walk and the Kirchhoff equations [Lj| allows to 
associate to each directed graph a networks of resistors 
(see again the Appendix) . In this framework, it is possi- 
ble to compute analytically the exponent a describing the 
growth of the resistance f2 as a function of the distance r 
between contacts, i.e. 51 ~ r a . Exploiting Eq. (TT2"j) . one 
can prove also the following Einstein relation 



a = 



2df_ 

d s 



(14) 



For local quantities and exponents, we can thus use 
Eq. (Ti"4"]) to calculate d s . We will verify it numerically, 
also testing the validity of the scaling hypothesis. 



B. Local spectral dimension 

When scaling holds, the local spectral dimension can 
be obtained using the scaling law of the resistance Q and 
the relation ([14")) . In the network of resistors associated 
to the Cantor graph every bidirectional link has resis- 
tance 1, while for unidirectional links, a whole bubble 
has unit resistance (see Appendix). The resistance be- 
tween the extremes of the structure of generation G can 
be evaluated using the recurrence relation: 

f2(G, n r ) = n r Q(G — 1, n r ) + (n r — 1). 

Using the initial condition f2(l,n r ) = 2n r — 1, one has 

0(G,n r ) = 2n°-l. (15) 

Expressing as a function of the distance r between con- 
tacts at generation G, in the G — > 00 asymptotic limit, 
we obtain the exponent a and, via the scaling relation 
(TT4")) . the value of d s , ruling the return probability and 
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the mean square displacement: 



0.5 n , , , , p 

n =2 n =3 d =1.22629 /■ 



i+ 



if n u > n r 



if T~L^ " T~by » 



if n u < n r 



r u s 



(16) 



and 7 = d s , where we have used that df = 1. Interest- 
ingly, the system exhibits a transition from a superdif- 
fusive to a normally diffusive regime at n u — n r . The 
random walker on the slim Cantor graphs experiences 
a superdiffusive regime while it presents a normal dif- 
fusion on fat structures, with logarithmic corrections in 
the critical case n u = n r . Intuitively, slim Cantor graphs 
mainly consist of long unidirectional bubbles, which lead 
to a superdiffusive behavior, while in the opposite case 
the bubbles practically disappear and hence one could 
expect a normal diffusion. 



C. Numerical results: The master equation and 
generalized scaling 

For a simple and effective numerical study of the sys- 
tem we solved iteratively the master equation Eq. @ 
with initial condition P^(0) = <%/2 and periodic bound- 
ary conditions. 

We first tested the dynamical scaling of the probabil- 
ity density (|12[) . considering for g(-) a periodic function 
with unit period in log„ (l(t)), so that the log-periodicity 
takes into account the self similarity of the underlying 
structure. It is thus convenient to look at the data at 
constant g(log„ u i(t)), namely, as prescribed by Eqs. (fTo) 
and (1161). at successive times 



h = t [r 



(17) 



Figs. 2] refers to a slim and a fat structure. For the 
former case the scaling is reasonably accurate, although 
the convergence to the asymptotic shape is not complete. 
Moreover the two panels show that for different sequences 
of times satisfying relation (fT?)) (i.e. different to) the scal- 
ing function changes, evidencing that the introduction of 
a generalized scaling with a log-periodic term is necessary 
in this situation. The scaling function also depends on 
the starting site i, while the exponent are site indepen- 
dent. For the fat case, the curves tend to be closer and 
closer to a Gaussian, confirming that the diffusion is nor- 
mal. Times have not to be chosen according to Eq. (TT7)) 
since in this case log-periodic oscillations are not present. 

Figs. [5] and [6] demonstrate that the probability of re- 
turn Pu{t) and the mean square displacement behave 
asymptotically (up to log-periodic corrections) as pre- 
scrived by Eqs. ©, JTOJ and (JTSJ (d/=l for all Can- 
tor graphs). Indeed, power-law fits of the curves give 
exponents 1.21, 1.13, 1.03 which are in excellent agree- 
ment with the theoretical values of d s , Eq. (ITBK ii, = 
1.2262 . . . , 1.11577 . . . , 1 respectively). 
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FIG. 4: Dynamical scaling of the probability for initial site 
i = 1. The size of the lattices correspond to generations 
G = 9, 8 respectively. 



IV. DIFFUSION ON CANTOR GRAPHS: 
AVERAGE BEHAVIOR 



The exponents characterizing the asymptotic behav- 
ior of random walks are in general independent of the 
starting site i, even on inhomogeneous structures . How- 
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FIG. 5: Dynamical scaling of the probability for initial site i = 
1. Lattices of size Ntot correspond to G = 9, 7, 8 respectively. 
To avoid the big jumps and zeros what is actually plotted is 
a coarse grained probability "Y^tLt Pu(t)/[tk — ifc-i]- 
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FIG. 6: Growth of the mean square displacements for initial 
site i = l. Lattices of size N corresponding to G = 9,8,7 
respectively. 



ever, due to inhomogeneity, averages over the starting 
points can result in radically different asymptotic behav- 
iors with respect to local case It is thus necessary 
to introduce some further averaging on the initial con- 
ditions. Let Sk (k integer) be a sequence of subgraphs 
(Sk C Sfc+i) covering the infinite graph. We define the 
average mean square displacement: 



1 



(s»)= hm ir>M) 



ies k 



(18) 



(TVfc is the number of sites in Sk) and the average return 
probability 



(19) 



Interestingly, due to the fact that the thermodynamic 
limit does not commute with the large time limit, the 
average quantities (|18|19|) can feature a different asymp- 
totic behavior with respect to their local counterparts. 
Therefore, as done for the local case, Eqs. ([SJ and (fT0|) . we 
introduce the exponents d s and 7 for the average quan- 
tities: 



(20) 



On the Cantor graph, average exponents not only differ 
from the local ones, but also they do not satisfy relations 
(JTTJ) and (fT4| . Indeed, for average quantities scaling is 
violated and the topology of the graph influences diffu- 
sion in a highly non trivial way. Diffusion is ballistic on 
slim (n u > n r ) Cantor graphs, while on fat graphs, even 
if d a always equals one, there are two different scenarios: 



7 = 



1 



log Tb u " ' U 

if < n r 



(21) 



In both cases, we obtain a lower bound for the average 
mean square displacement and then we check numeri- 
cally that this lower bound is satisfied when n r < n^, 
obtaining the exponent given in Eq. (|21[) . 

In analogy with [8|, let us study the probability Pb(l) 
to take a first ballistic step of length > I. We will show 
that, on slim Cantor graphs, this probability approaches 
1 when G — > 00 for every I, and the motion is therefore 
ballistic on average. On the other hand, on fat Cantor 
graphs, the limit is a finite number going to for I — > 00 
with a characteristic power law. Let us consider the uni- 
directional segments of length greater than I: the "favor- 
able" sites in a unidirectional segment of length Lk > I 

are Lk — l, and we have to sum over k. Recalling that the 

,G-i-fc 



number of segments of length Lk 
we can write Pb{l) as 



i£ is (n r - 1)? 



Ett, K - 1 - K - 1) K) 

N tot (G 7 n s ,n u ) 



G-l-k 



P b {l)= hm 

where k is an integer chosen such that > I, i.e. k := 

(reminding that \x~\ := min{n G Z | n > x} 

is the ceiling function), and N tot (G,n r ,n u ) is the total 
number of sites in the graph of generation G, given in 
([7]). The series in the numerator of (|2"2"|) can be summed 
and we obtain: 



log(0 
log n„ 



Pb(l) 



1 



lim . 

G->oo N u 



n r — 1 



+ 



n r - 






n r — 


n u 


\n r ) 



, — ko 



nf + (l + l) 



ies k 



Therefore, we have two exponential terms with different 
bases: the leading term for G — >• 00 will then be different 
depending on n r and n u . 
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If n u > n r , ( i.e. on a slim Cantor graph), considering 
the asymptotic behavior of N tot (G, n r ,n u ) in Eq. (JTJ), one 
obtains 



Pb(D 



<K-i) + o(<) 

<££o - l) + o(nG) 



lim 



1. 



Thus, the probability of going through a ballistic step 
of length > I is 1 for every / and the walker behaves 
ballistically on average. 

Let us now consider the same limit in the case n u < n r . 
Omitting the ceiling function in kg (it can be done with 
a suitably chosen Z), we have 



n(0 



—fen 



z 1 



(23) 



In order to estimate the average mean square displace- 
ment in the fat case, let us introduce the average proba- 
bility of being at distance r after t steps as P(r,t). It is 
convenient to split it in two parts writing 



P(r,t)=P*(r,t;r < t) + 



Pb(t) 



(S(r-t) + S(r + t)) 



where P{> is the average probability of performing t con- 
secutive ballistic steps and P*(r,t;r < t) is the average 
probability of arriving in t step at r after some scatter- 
ing. Eq. ([24| provides a lower-bound for the mean square 
displacement 



dr r 2 P(r, t) 



> 



(25) 



where we used expressions ([23)) for P&(£). For 1 < 
^ < 2 the inequality (|25|) proves that the system is 

o_ log n r 

superdiffusive and (x 2 ) ~ t lo s»" is expected to be the 
correct asymptotic behavior of the average mean square 
displacement, In Fig.[7]we compare our analytical predic- 
tion with numerical data, showing an excellent agreement 
[2ll ] . For ™ r > 2 inequality (|25|) is trivial and a normal 
diffusion is expected, since the ballis tic st retch does not 
provide a significant contribution to (ir 2 ). The presence 
of normal diffusion in this case is also evidenced in Fig. [7] 
Fig. [8] shows that, in the fat case, there are actually 
two different contributions to P(r,t): a ballistic peak 
lowering with t and a central peak which scales as in 
normal diffusion. This behavior clearly breaks the scal- 
ing hypothesis (|12l) and relations (ITT1) and (ITi|) between 
exponents do not hold. Indeed, for log(n r )/ log(n«) < 2, 
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log rii 



, while d s = 1 and a = 1. 



CONCLUSIONS 



The recent experiments on light scattering in disor- 
dered materials open the way to direct and tunable 




FIG. 7: Average mean square displacement, as denned in 
Eq. (|18[) f° r the fat fractals with n r = 4, n u = 3 and 
n r = 5, n„ = 2. Solid and dashed lines are the superdif- 
fusive, Eq. (I25p . and diffusive behaviour expected in the two 
structures, respectively. 
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FIG. 8: Dynamic scaling of the averaged probability P with 
n r = 4 and n„ = 3 lattices of size N to t = 58976 (G — 7). The 
average is on all initial sites. 



measurements on systems performing Levy walks. The 
quenched disorder in the experimental samples has been 
engineered in order to obtain a specific effects, and the 
geometry induces a correlation in Levy walks. In this pa- 
per, we take a first step in understanding the influence of 
geometry in these disordered samples by studying a one 
dimensional structure obtained from the generalization of 
Cantor and Cantor-Smith- Volterra sets, where a random 
walker performs steps that are both Levy distributed 
and correlated. Although the applicability of such ide- 
alized one-dimensional models to real experiments is not 
granted, much theoretical insight can be obtained from 
their solution. Using an exact mapping to an equivalent 
directed graph for the process, we determine the asymp- 
totic behavior for the return probability and the mean 
square displacement and study the validity of the Ein- 
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stein relations in the local and average case. 

From our results, it appears that the average case is 
the closest to the real experimental setting 0, and fur- 
ther analysis in higher dimensions are currently under 
consideration. An important point is that experimental 
samples do contain a certain degree of disorder, which 
must be taken into account in order to understand to 
what extent it can affect diffusion properties. 

Appendix 

In this Appendix we show that the Einstein relation 
Eq. (fT4f between different exponents is a direct conse- 
quence of the generalized scaling hypothesis, Eq. (fi~2")) . 

Since Pij(t) is a probability, its sum over all sites n 
is normalized to one at every time t. Moreover, in the 
scaling hypothesis, all the spatial dependence is encoded 
by the distance Xij = r, therefore such a sum can be 
evaluated by introducing an integral in Kr df ~ 1 dr where 
df is the fractal dimension and AT is a suitable constant, 
obtaining 

t~ d °l 2 j h(^,g{\ogJ{t))\Kr d f- l dr = 1. (A.I) 
Changing the integration variable to r/£{t) we have 

t -d s /2 fdf^ G ^ = j ( A 2 ) 

where G(t) is a log-periodic function. Eq. (|A.2[) is satis- 
fied if G(t) = 1 and if Eq. (T5J} holds. 

To derive Eq. (fTT|) . let us consider the mean square 
displacement: 

<*<(*)> =t~ ds/2 J r 2 U (^,g(log x £(t))j Kr^-Ur 

(A.3) 

Upon choosing again the new integration variable r/£(t), 
we have 

(x 2 ) = £{t) 2 Gi(t) = &' d i G x (t) (A.4) 

where G\ (t) is a log-periodic function. 

Let us finally consider the equation for the electric po- 
tential V; on a network of unitary resistors where a uni- 
tary current flows from site to site n, we have: 

-^LjiVj =5 i0 -S in (A.5) 

3 



where Lij = ZiSij — Aij is the Laplacian matrix. Notice 
that according to (|A.5[) . in a directed Cantor graph the 
inner links of a bubble are short-circuit, and a whole 
bubble has resistance one. In the framework of random 
walks on directed Cantor Graph the master equation ([2]) 
can be recasted as: 

P 0i (t + 1) - P oi (t) = -J2 L j,i p oj{t)/zj + S l0 S t o- (A.6) 

3 

Denoting with Poi{uj) the Fourier transform of Poi(t) we 
get 

P 0i (u)(e iw - 1) = - J2 L i,iPoi(u)/zj + 5 i0 (A.7) 

3 

comparing Eqs. ([A~5]) and (|A~7| 

Vi = - Km (Poi(w) - P ni (u)) (A.8) 

The potential difference between sites and n as a func- 
tion of their distance r, can be, hence, obtained intro- 
ducing in lim^-i.0 Pon{u>) the scaling relation Eq. (|12p 

V{r)~hmJ e^r^/Vif^.sOo&A*))) (A.9) 
Changing the variable of integration into t' — cot we get 
V (r) ~ r 2d f^- d f Hm G 3 (j-, g(log x £(cj))) (A.10) 

where G3 is a suitable function and £(u>) = Lo~ ds ^ 2d ^ is 
the correlation length in terms of the frequency uj. There- 
fore, Eq. (JHJ) holds. 
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